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ABSTRACT 

We identified 46 unresolved source candidates in the Hubble Ultra Deep Field, 
down to 1775 = 29.5. Unresolved objects were identified using a parameter S, 
which measures the deviation from the curve-of-growth of a point source. Exten- 
sive testing of this parameter was carried out, including the effects of decreasing 
signal-to-noise and of the apparent motions of stars, which demonstrated that 
stars brighter than 1775 = 27.0 could be robustly identified. Low resolution grism 
spectra of the 28 objects brighter than i 775 = 27.0 identify 18 M and later stellar 
type dwarfs, 2 candidate L-dwarfs, 2 QSOs, and 4 white dwarfs. Using the ob- 
served population of dwarfs with spectral type M4 or later, we derive a Galactic 
disk scale height of 400 ±100 pc for M and L stars. The local white dwarf density 
is computed to be as high as (1.1 ±0.3) x 10~ 2 stars/pc 3 . Based on observations 
taken 73 days apart, we determined that no object in the field has a proper mo- 
tion larger than 0.027" /year (3a detection limit). No high velocity white dwarfs 
were identified in the HUDF, and all four candidates appear more likely to be 
part of the Galactic thick disk. The lack of detected halo white dwarfs implies 
that, if the dark matter halo is 12 Gyr old, white dwarfs account for less than 
10% of the dark matter halo mass. 
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1. Introduction 

Studying the stellar content of our galaxy down to the faintest possible magnitudes is 
necessary to study the structures of the Galactic disk and population II halo (referred to 
as Galactic halo thereafter) and stellar population (Mendez et al. 1996). The nature of the 
Galactic dark matter halo has in addition been the subject of some debate since the MACHO 
project (Alcock et al. 2000) observed lensing events in the direction of the Large Magellanic 
Cloud (LMC) and claimed that 20% of the dark matter halo mass could be in the form of 
low mass objects with masses smaller than l.OM (Alcock et al. 2000). While there are 
alternative explanations for the observed lensing events (Sahu 1994; Sahu 2003; Gates & 
Gyuk 2001, and references therein), the detection of faint, blue unresolved objects with high 
proper motion has been used as evidence for the existence of a large population of old white 
dwarfs in the Galactic dark matter halo. Oppenheimer et al. (2001) estimated that about 
2% of the dark matter halo mass could be accounted for by old halo white dwarfs. There 
has since been some debate as to whether the former result might have been caused by a 
dynamically warmer thick disk component (Reid, Sahu & Hawley 2001) and the latter might 
have been affected by the misidentification of faint blue extra-galactic sources (Kilic et al. 
(2004), also see Reid (2004) for a complete review). The existence of ultra cool white dwarfs 
(T ^ 4000° K) has however been confirmed in recent SLOAN data (Gates et al. 2004). 
A search for distant white dwarfs using deeper observations provides another opportunity to 
firmly examine the existence of a significant population of old white dwarfs in the Galactic 
halo and the dark matter halo. As shown by Kawaler (1996), deep broad band imaging can 
probe a significant portion of the Galactic dark matter halo and should identify Galactic 
dark halo white dwarfs in the process. Ideally, one would like to identify white dwarf candi- 
dates spectroscopically, instead of relying on broadband colors only, and to determine their 
Galactic disk or halo memberships using sensitive proper motion measurements. 
The Advanced Camera for Survey (ACS) HST observations of the Hubble Ultra Deep Field 
(HUDF, Beckwith et al. 2004) and their unprecedented depths offer a new opportunity to de- 
tect faint stars, such as white dwarfs, at high galactic latitude. Additionally, the GRAPES 
survey (GRism ACS Program for Extragalactic Science. PI: Malhotra, see description in 
Pirzkal et al. 2004) provides low resolution spectra of most sources in the HUDF. These 
spectra can be used to confidently determine spectral types of stars in the HUDF, down to 
i775 = 27.0 mag. More importantly, these spectra allow one to differentiate unresolved, blue, 
faint extra-galactic objects (e.g. QSO) from white dwarfs and other stars. 
Section 3 outlines how unresolved sources were identified in the HUDF, while the spectro- 
scopic identification of unresolved sources is discussed in Section 4. Section 5 describes an 
analysis of the cumulative number distribution of the late-type dwarfs identified in the HUDF 
and its implication for the Galactic disk structure (Section 5.1). Section 5 also includes a 
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description of the white dwarf content of the Galactic disk and halo (Section 5.2). 

2. Description of the data 

The HUDF was imaged using the ACS and the F435W (B 435 ), F606W (V 606 ), F775W 
(i 775 ) and F850LP (z 850 ) bands at different epochs. The observations were carried out by 
cycling through each filter while using a small dither pattern to mask out the ACS inter-chip 
gap and other cosmetic problems. The HUDF was observed during 4 different visits, rotating 
the field of view to a different position angle each time (0°, 8°, 85°, 91° with respect to the 
first visit, Beckwith et al. 2004). The individual ACS images had an integration time of 
1200 seconds, a size of 4096x4096 pixels and a pixel scale of 0.050" /pixel. These images were 
combined together by drizzling them onto a finer pixel grid with a pixel scale of 0.030" /pixel, 
while geometrically correcting these images. The limiting AB magnitudes (S/N=10) of the 
i 775 and z 850 images are 29.0 and 28.7 mag, respectively. In this paper, images taken from 
09/24/2003 to 10/28/2003 and from 12/04/2003 to 01/15/2004 are referred to as Epoch 1 
and 2 images, respectively. The elapsed time between the mean observing dates of the two 
epochs is 73 days. Images from Epoch 1 and Epoch 2 were taken at position angles that are 
nearly 90° apart. 

While the deep, combined ACS HUDF images described above were used to determine which 
objects are unresolved (Section 3), eight special partial image stacks were also assembled to 
take advantage of the fact that the data were taken over a period of several months. Four 
partial image stacks were created for each HUDF Epoch, using the exact same method used 
to create the original combined HUDF images. These partial image stacks can be used to 
determine the variability (Cohen et al. 2004) as well as measure the proper and parallax 
motions of the objects in the field. The eight partial image stacks have the same physical 
properties as the combined ACS HUDF images but with one eighth of the integration time. 
Note that the cosmic ray rejection step was done separately while assembling each partial 
image stack so that any strongly varying or moving object would not be removed by the 
cosmic ray rejection process. The astrometric registration of the combined HUDF images 
and the eight partial image stacks is estimated to be better than 0.1 ACS pixel (0.005" or 
0.16 HUDF image pixel) (Beckwith et al. 2004). 

Spectroscopic observations of the UDF using the ACS slitless grism mode (G800L) were 
carried out as part of the GRAPES project. The reader is referred to Pirzkal et al. (2004) for 
an in-depth description of these spectra. These spectra are low resolution spectra (R=100), 
span a wavelength range of 5500A ^ A ^ 10500A, and can be used to spectroscopically 
identify objects in the field. 

The HUDF Public Catalog 1.0, containing 10181 objects, was used to identify objects in the 
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field and all objects are referred to by their ID in this catalog (UID). The magnitudes of 
these objects were re-measured using matched aperture photometry with apertures defined 
using Sextractor (Bertin & Arnouts 1996) and the i 775 band HUDF image. The photometric 
ACS magnitude zero-points (25.65288, 26.49341, 25.64053, 24.84315 in the B 435 , V 606 , i 77 5, 
z 8 5o bands respectively) were taken from the GOODS survey (Giavalisco et al. 2004) . 



3. Identification of point sources in the HUDF 

We identified unresolved objects in the HUDF using a single-index (S) method that is 
based on a simple analysis of the curve of growth of the light distribution of an object. The 
IRAF task RADPROF was used to produce these light curves of growth. This task first 
computes the intensity weighted center of the object (in a 3x3 pixel box), and then performs 
aperture photometry using increasing aperture sizes (up to a radius of 10 pixels). These 
photometric measurements are then fit (using eight 3 a rejection cycles) to a cubic spline 
of order 5. RADPROF was applied to stamp images of objects which were masked using 
Sextractor segmentation maps. These segmentation maps were used to avoid contamination 
from neighboring objects. The RADPROF parameters described here were selected because, 
following extensive initial testing, they were found to result in curves of growth that are not 
too sensitive to the presence of faint nearby objects while consistently remaining sensitive to 
the presence of faint, diffuse extended emission around objects. The stellarity S of an object 
is defined as 

S 2 = ~ £(i^(r) - F psf {r)f 

i rmax r=Q 

where F Q bj(r) and F ps f(r) denote the curve of growth of light, normalized to 1.0 at a distance 
of r max , of the object and of a point source. A value of 10. was used for r max . 

The index S is then simply the square root of the sum of the squared differences between 
the two curves. Unresolved objects have S values that are close to 0.0 while increasingly 
resolved objects produce increasingly larger (unbound) values of S. The detailed calibration 
of S was done in several stages, each time refining the reference curve of growth that was 
used. This is described in some detail below. 

First, a few obvious, bright, unresolved and unsaturated sources were pre-selected for further 
study. This was done in both the and z 8 5o band images separately. The GRAPES spectra 
of these sources were examined and a main-sequence star was selected (UID 911). The curves 
of growth of this object were measured separately in the i 775 and z 850 bands and were used 
to compute S\ 7n and S Z850 values for all 10181 objects in the HUDF. A new set of 25 bright 
point sources, well distributed over the whole ACS field, were then carefully selected based on 
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their low S i775 values (arbitrarily chosen to be 5 zg50 < 0.04). The GRAPES spectra of these 
25 objects were visually inspected to ensure that these objects were stars. Stamp images of 
these objects were individually drizzled (Fruchter & Hook 2002) and carefully re-centered 
to produce a stack of images which were then combined together to produce empirical 1775 
and z§5o PSF images. A sigma-based rejection algorithm was used to avoid the effect of 
any faint neighbor in the final combined images. Final 1775 and zgso band reference curves 
of growth were generated using these empirical PSF estimates and all of the S i775 and <S Z850 
values were re-computed. Figure 1 shows the distributions of the S values in both the 1775 
and z 8 5o bands, as a function of magnitude. There are clearly two distinct groups of object 
in the HUDF images: one group consists of object that are completely or partially resolved, 
forming a large distribution around the S value of 0.2 and 1775 = 28.0 mag. A second group 
consists of unresolved objects with low S values at all magnitudes. Bright stars in Figure 1 
are observed to have S values which are consistently lower than 0.05. The exception to this 
are bright stars near i7 75 = 18 mag. which are saturated and have artificially high values of 
S (~ 0.1) . The effect of bright star saturation on the measured values of S were investigated 
and values of 0.1 ^ S < 0.15 are expected if the center of a star is saturated or missing in 
the HUDF images. Nevertheless, the populations of resolved and unresolved objects appear 
to be well separated at brighter magnitudes in Figure 1. 

Unresolved sources in the HUDF were selected using the criterion that <Si 775 ^ 0.05 mag. Four 
bright saturated stars were added to the final object list manually. The final list of unresolved 
UDF sources contains 46 objects, spanning the magnitude ranges of 19.0 ^ i 775 ^ 29.5 mag 
and 18.3 ^ z 85 o ^ 29.2 mag. Table 1 summarizes some of the properties of these objects. 

3.1. Testing the selection process 

The ability to distinguish between resolved and unresolved objects using S values was 
extensively tested, using a subset of the HUDF object catalog containing the 369 brightest 
( z 850 ^ 25.0 mag) objects in the field. Simulated 1775 and z 850 band images of these objects 
- made 0.0, 0.75, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, and 9.0 magnitudes fainter — were created 
and S values measured using each simulated image, exactly as described in Section 3. The 
S value of each of these 369 test objects was thus measured 10 times, each measurement 
corresponding to a different apparent magnitude. Figure 2 shows how the S value remains 
largely unchanged as an object is made fainter and noisier, all the way down to i 775 = 29.0 
mag. The vertical pattern seen in this figure is a result of the 10 magnitude steps used in this 
simulation. This test demonstrates a very small fraction of misclassification: not a single 
extended object was misclassified as a star (e.g. no object with S ^ 0.05 originally in the 
HUDF images was measured to have S ^ 0.05 as the object was made fainter), and fewer 
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than 5% of stars were misclassified down to i 775 = 29.0 mag. While some objects shown in 
Table 1 are fainter than this, only objects significantly brighter than this limit (1775 ^ 27.0 
mag) are actually included in the analysis described below and these are expected to be bona 
fide unresolved objects. 

3.2. Testing the effect of proper motion 

Some of the stars in the HUDF objects could have a high tangential-velocity (V t ) and a 
correspondingly large proper motion. Such objects would be blurred in the combined HUDF 
images and would appear to be resolved. Similarly, nearby objects would suffer from a 
significant amount of parallax, again causing these objects to be blurred and mis-identified. 
The index S proved relatively insensitive to a small amount of object motion. This is not 
unexpected since S was tuned to the detection of faint extended components, The curves of 
growth produced by RADPROF were computed by summing, using several sigma rejection 
steps, the pixel intensities inside annuli of increasing radii. The effect of a nearby source 
is strongly diminished as its separation from the main source increases and the number 
of unaffected pixels increases. In addition, the use of a fifth order spline imposes some 
smoothness on the curve of growth which further lowers the effect of a nearby object. Still, 
the effect of object motion on measured values of S was carefully investigated by generating 
images of the 46 unresolved HUDF objects that were selected as described in Section 3. 
These were appropriately shift-and-added, accounting for the orientation and pointing used 
for each individual HUDF observation while allowing objects to move as a function of time, 
to simulate the effect of proper motion. As larger values of proper motion were introduced, 
the S values increased only slightly and S was found to be mostly unaffected by shifts smaller 
than 2.5 HUDF pixel (fx = 0.3 "/year). 70% of objects were still properly identified as point 
sources when shifts 2.5 and 5.0 HUDF pixels were applied (0.3 ^ \i ^ 0.6). 

4. Identifying stars 

4.1. Using Grism spectra to identify stars 

The spectral type of 28 unresolved objects brighter than i775 = 27.0 could be determined 
robustly using template fitting of their GRAPES spectra (Pirzkal et al. 2004). The template 
set included Pickles (1998) main-sequence templates, L and T dwarf spectra (Kirkpatrick 
ct al. 2000; Burgasser et al. 2003), DA and DB white dwarf spectra (Harris et al. 2003) 
obtained from the SDSS archive, and some low temperature (3000°K, 5000°K, and 8000°K) 
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model spectra of white dwarfs (Harris et al. 2001). The quality of each spectral fit was 
evaluated visually and the results of these fits are given in Table 2. Objects too faint to have 
reliable fits are indicated. The GRAPES spectra and the two best-fits listed in Table 2 are 
shown in Figures 3, 4, and 5. Unresolved extra-galactic sources are labelled as GAL in Table 
2 and are discussed in more details in Section 4.3. 

Of the 46 point sources detected in the HUDF, 26 stars, all brighter than i 775 = 27.0 mag 
were identified. As shown in Table 2, the majority of stars in the HUDF (18 out of 26 ) 
are found to be M dwarfs. Figures 6 and 7 show these objects in B 435 vs. B 435 -V 6 o6 and 
i775vs. i775-V606 plots. Such plots have in the past shown a lack of faint red objects in the 
region V 60 6 ^ 27.0 mag and (B 435 — V 60 6) ^ 1-0 (Mendez et al. 1996). Objects found in this 
region of the plot, appear to be either extra-galactic or L-type dwarf spectra (UID 443, 366). 
Further NIR observations will be required to verify the nature of these 2 L-dwarf candidates 
by confirming that they both have the expected blue NIR colors of this type of object. 
The HUDF field is at high galactic latitude (b=223°, l=-54°) and the effects of extinction 
and reddening are small (Ay = 0.026, Schlegel et al. 1998). The distances to the HUDF stars, 
which are shown in Table 2, were determined using their observed V 60 6 magnitudes, their 
estimated stellar types, and the dwarf star main-sequence absolute Johnson V magnitudes 
from Allen (2000). The latter were converted to the V 6 o6 filter bandpass and AB magnitudes 
to be consistent with the data at hand. Distance estimates obtained under the assumption 
that these objects are white dwarfs with intrinsic magnitudes of M mv = 14.0 and 17.0 are 
also given. 

4.2. Proper motion measurements 

The Galactic halo is isotropic and has no or little overall rotation. Halo objects are 
therefore expected to lag behind the rapidly rotating Galactic disk component (and the local 
standard of rest) with an average velocity of 200 ± lOOkm/s . As done in Oppenheimer 
et al. (2001), it is assumed here that Galactic halo objects have a velocity that is more than 
2o" above the velocity expected for an object that is in the Galactic disk (i.e. ^ 100 km/s, 
Chiba et al. 2000). 

The HUDF image stacks (Section 2) were subtracted from one another and the results were 
examined to see if any object, whether classified to be extended or not based on its <S Z850 
value, had moved significantly (e.g. more than 2.5 UDF pixel) from one image to the next. 
No such object was found. 

The motion of objects in the HUDF was further investigated by carefully measuring the 
position of objects separately in each of the eight available i775 band HUDF partial image 
stacks. These positions were measured using the IDL task MPFIT2DPEAK 7 , using a fixed 
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width (FWHM=2.5 pixel) Gaussian fit. Initially, the positions of the 360 most compact 
sources (S ^0.1 and i 775 ^ 28.5 mag) were examined. An average position and a standard 
deviation of the mean were computed for each object in the field using Epoch 1 and Epoch 
2 HUDF image stacks separately. Average Epoch 1 and Epoch 2 positions were subtracted 
from one another to compute the observed shift of each object. These are shown (scaled up 
by a factor of 500) in the left panel of Figure 8. Repeated measurements of the position 
of sources in the field using HUDF image stacks taken within the same Epoch were always 
within 0.1 HUDF pixel, or 0.003". Still, and as shown in Figure 8, there is a systematic, 
field dependent disagreement between the Epoch 1 and Epoch 2 measurements. This is likely 
caused by Epoch 1 and Epoch 2 images having been taken at position angles that are 90° 
apart, and/or the effect of telescope breathing, and/or very small amount of residual image 
distortion that is not taken into account in the current model of the ACS distortion maps. 
To produce a 3a detection with these data, an object must have a measured proper motion 
larger than 0.54 HUDF pixel or fi ^ 0.081 "/year. No object in Table 1 was observed to 
have such a large proper motion. 

The situation was however improved using a simple third order two dimensional polynomial 
fit to the observed distribution shown in the left panel of Figure 8. This fit, when applied as 
a correction, results in significantly lower systematics across the field, as shown in the right 
panel of Figure 8. The accuracy of the corrected shifts (of compact objects) between Epoch 
1 and Epoch 2 is improved by a factor of 3, with a 3a detection now corresponding to a shift 
of 0.18 HUDF pixel or /i > 0.027 "/year. 

The search for high proper motion objects in the HUDF was extended to all 10181 objects in 
the field by measuring and correcting (using the correction determined above) the observed 
shifts of these objects. Only four sources were found to have moved significantly (but with 
a significance slightly lower than 3a) from Epoch 1 to Epoch 2. These objects are UID 443, 
9020, 911, and 7525 and they are all unresolved objects that are listed in Table 1. Their 
measured shifts are 0.17, 0.17, 0.14, and 0.15 HUDF pixels respectively, or approximately 
0.02" /year, which is smaller than our imposed 3a limit. 

4.3. Extragalactic Unresolved Sources 

Two unresolved, extra-galactic sources with i 775 ^ 27.0 were identified (UID 6732, 9397). 
The GRAPES spectrum of these objects show prominent emission lines and the redshift of 
these objects is estimated to be z=3.2 and z=3.0, respectively (Xu 2004). Identification of 
such objects using GRAPES spectra could easily be carried out down to magnitudes fainter 



7 http://astrog. physics. wisc.edu/ craigm/idl/ 
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than i 77 5=27.0 mag, as demonstrated by the fact that several fainter extragalactic objects 
with 27.0 < i 775 < 29.5 were easily identified (e.g. UID 4120, 8157). Without spectral 
information, these objects would have been misidentified as stars (e.g. UID 9397 would be a 
good candidate for a cool white dwarf at a distance of ~ 1 Kpc) . The issue of contamination 
by faint, blue extra-galactic objects in previous studies (Ibata et al. 1999; Mendez & Minniti 
2000) has been discussed in the past (Kilic et al. 2004) and Figures 6 and 7 illustrate how 
these objects could be mistaken for a variety of stellar objects ranging from M and L dwarfs 
to white dwarfs. Spectroscopic confirmation is therefore crucial to accurately determine 
stellar counts in deep fields such as the HUDF. There are 2 extra-galactic objects amongst 
the 28 point sources with i 775 ^ 27.0 mag, which corresponds to a 7% contamination level 
at these Galactic coordinates. 

The 2 extra-galactic objects (UID 6732, 9397) were used to check the accuracy of the proper 
motion measurements described in Section 3.2. The measured shifts of these objects between 
Epoch 1 and Epoch 2 are smaller than 0.04 HUDF pixel (0.0012"), which is what one would 
expect based on the proper motion error estimate quoted in Section 3.2). 
There are a total of 14 faint (27.0 < i 775 < 29.5), blue ((B 4 35-V 6 o 6 ) ^ 1.8, (V 606 -i 7 75) < 0.8) 
objects listed in Table 1 that were not included in the analysis above because they are either 
too faint to have useable GRAPES spectra or robust S based classifications. Some of these 
sources might well be faint white dwarfs in either the Galactic disk or halo but are hard 
to distinguish from extra-galactic sources without higher signal-to-noise spectra or more 
accurate proper motion measurements. 

5. Discussion 

5.1. M and L Dwarfs and the Galactic disk scale height 

18 objects were identified to be M (or early L) dwarf main-sequence stars. Two objects 
were found to be best fit by early L dwarf templates (Kirkpatrick et al. 2000). Both of these 
objects are red enough to be outside of the color range shown in Figures 9 and 10 with 
(V606 ~ 1775) = 2.21 and (V606 _ 1775) = 3.47. Deep infrared observations of these two objects 
will be required to confirm that they are L dwarfs. NICMOS observations of the HUDF did 
not detect these objects in either the J or H bands. All other very red objects are best fit 
using M dwarf templates. The broad spectral features of M dwarfs listed in Table 2 were 
well fitted using M dwarf Pickles (1998) templates, and for objects brighter than i 775 = 27.0 
mag we believe the identification to be secure. 

Because these late-type sources can be seen to great distances, they can be used to probe 
the structure of the Galactic disk and halo as traced by the lowest-mass stellar components 
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of these populations. This was explored using Monte Carlo mass function (MF) simulations 
based on those developed by Burgasser et al. (2004). Assuming power-law representations 
of the MF, dN/dM oc M~ a , for masses 0.005 to 0.2 M , luminosity (LF) and effective 
temperature {T e ff) distributions were created using evolutionary models from Burrows et 
al. (1997) over the range 100 < T e // < 4000 K. For all simulations, a was fixed at 1.13 for 
0.1 < M < 0.2 M (Reid et al. 1999) but allowed to vary between 0.5-1.5 for M < 0.1 M ; 
number densities are normalized to the empirical value at 0.1 M from the 8 pc sample (Reid 
et al. 1999). 

Two separate populations were simulated: a disk population with a flat age distribution 
spanning 10 Myr to 10 Gyr (the majority of sources have age greater than 1 Gyr), and 
a halo population with uniformly sampled ages between 9-10 Gyr and a relative number 
density of 0.25% (Digby et al. 2003). The resulting bolometric LFs were converted to / 
(Cousins) LF's using a polynomial fit to empirical data (spanning 12.4 < Mi < 22.2; i.e., M4 
to T8) from Dahn et al. (2002) and Golimowski et al. (2004). Note that this fit is based on 
measurements for Solar metallicity field dwarfs, and may not be appropriate for a subsolar 
metallicity halo subdwarf population. Apparent /-band distributions were determined by 
first assuming a constant density population out to a limiting magnitude m/ = 27.5 (taking 
into account the difference between HUDF AB-magnitudes and Vega magnitudes from the 
empirical data) and then applying a correction for the vertical distribution of sources. For 
the disk population, we assumed a density distribution n cx sech(zfho) 2 (Reid & Hawley 
2000), where z = dsin/3, d is the source distance, (3 = 54.5° is the Galactic latitude of the 
HUDF, and h is the 1/e disk scale height, assumed to range between 200 and 500 pc. For 
the halo population, a Galactocentric oblate spheroid distribution as given in Digby et al. 

(2003) was assumed, using the values derived there and an axial ratio q = c/a = 0.7. 

Figure 11 presents the results of these simulations in the form of cumulative number 
distributions as a function of i 7 j 5 magnitude down to the limiting magnitude of our sample. 
These distributions show that variations in the disk scale height are far more pronounced 
than those from the different MF's assumed, and this analysis is limited to the former 
parameter. The cumulative distribution of the six dwarfs in Table 2 with spectral types 
M4 and later (consistent with the mass constraints of the simulations) is also shown. The 
observed distribution matches that of the h = 400 pc disk simulation very well, particularly 
out to i775 < 23.5 mag where the halo population makes negligible contribution. Because the 
simulated distributions are rather sensitive to the disk scale height, particularly at fainter 
apparent magnitudes, one can conservatively constrain h to ±100 pc assuming Poisson 
uncertainties. This is in good agreement with the estimate of 340 ± 84pc from Ryan et al. 

(2004) which was computed using ACS HUDF parallel fields i-drops. For i 775 ^ 25 mag, the 
observed distribution is slightly greater than the disk population alone but does not increase 
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as sharply as the combined disk+halo distribution. Indeed, the shape of the observed number 
distribution implies few if any halo stars in the HUDF down to i 775 = 27.0 mag, consistent 
with the lack of significant proper motion sources in this sample. This suggests that either 
the number density ratio of halo to disk stars, or the adopted axial ratio for the halo density 
distribution, or both, may be smaller than assumed here. With no late-type halo subdwarf 
detections in this sample, we cannot usefully constrain these possibilities. 

5.2. White Dwarfs 

5.2.1. White Dwarf candidates 

Objects brighter than i 775 = 27.0 were individually examined and distances were esti- 
mated under various assumptions of what the exact nature of each object might be: (1) a 
main sequence star; (2) a young white dwarf; (3) an old white dwarf. The reddest objects, 
which are all well-fitted by M-and-later type templates, as well as the few extra-galactic 
objects identified in Section 4.3 were excluded from this analysis. The remaining 8 objects 
were all initially considered to potentially be white dwarfs. 

Assuming that an object is a white dwarf, one can use the measured V606 _ 1775 color 
of that object (Table 1) and the cooling curves of Richer et al. (2000) (after accounting for 
passbands and zeropoints differences), to derive an absolute Vgo6 band magnitude, an age, 
and a distance to that object. As shown in figure 2 of Richer et al. (2000), blue white dwarfs 
can either be young, hot objects (Age ^ 5 Gyr, T 10 4 K, My 606 ^ 16.0 mag) or older, 
cooler, intrinsically dimmer ones (Age ^ 10 Gyr, Tr;3x 10 3 K, My 606 ^ 16.0 mag). This 
results in several white dwarf distance estimates for each object. 

Proceeding via a process of elimination, the distinction between main sequence stars, young 
white dwarfs, old white dwarfs, and disk or halo white dwarfs is possible. Assuming that a 
star is a main sequence star, which is intrinsically much more luminous than a white dwarf, 
can cause the distance estimate for that object to be unreasonably large. Assuming that a 
star is a young white dwarf (Age ^ 5.0 Gyr) implies that this object is less likely to be part 
of the Galactic halo since the latter is composed of much older objects. Also, and based on 
the lack of proper motion detection discussed in Section 4.2, one can define an upper limit to 
the tangential velocity (Vr m ox) °f eacn object once its distance is determined. The motion 
of an object moving by more than 0.18 HUDF pixel (0.027 "/ yr) should have been detected 
at the 3a significance level as discussed in Section 3.2. As discussed above, objects in the 
Galactic halo are expected to have velocities around 100-200 km/s. Even if projection effects 
should result in lower values of Vr, one would expect at least some of them to have Vr values 
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larger than 30-60 km/s. A final clue to help narrow down the nature of a particular object 
is provided by the fact that an object that is nearby would suffer from a significant, easily 
measurable, parallax during the 73 days interval between Epoch 1 and Epoch 2 observations. 
Calculating the parallax vector for the HUDF exposures shows that an object 30 pc away 
would produce a parallax of about 0.060" or 2 HUDF pixel between Epoch 1 and Epoch 
2 observations. Objects that are closer than 200 pc are therefore expected to produce a 
parallax that would be detected at more than 3<r level. 

Examining Table 3, objects 4322, 4839, 7768, 9020 are unlikely to be main-sequence 
stars as this would make them very distant objects. The Magellanic Stream (Mathewson 
& Ford 1984) is too far away from the HUDF to possibly explain the existence of main- 
sequence stars at such large distances. Objects 1147, 3166, and 5921 which do not (this 
is true of all objects listed in Table 3) have measurable parallax or proper motion, are not 
likely to be white dwarfs. Being either young or old white dwarfs would place these two 
objects less than 200 pc away and would have resulted in measurable parallax in the HUDF 
image stacks. Another object unlikely to be a white dwarf is object 9230 which was ob- 
served to have a fainter, unresolved companion 0.5" away. This companion, for which we 
have no GRAPES spectrum, has a V 60 6 band magnitude of V 60 6 = 23.7 mag and the colors 
B435 — V 60 6 = 1.6, V 60 6 — 1775 = 0.9, i 775 — z 8 5o = 0.4 mag. Based on these colors, this object 
should be an early M dwarf that is (4.9 ± 1.2) x 10 3 pc away. This places both object 9230 
and its companion at the same distance (within error bars). The likelihood of any two stars 
being within 0.5" in the HUDF, and at nearly the same distance, is very small (^ 5 x 10~ 3 ) 
and these two objects are likely to be part of a binary system, with object 9230 being a K 
type main-sequence star. Independently of this, object 9230 cannot be an old white dwarf 
since this would place this object at distance of 34 pc where its parallax motion would have 
been very easily detected. 



5.2.2. Disk or halo white dwarfs? 

Four objects (4322, 4839, 7768, 9020) remain as white dwarf candidates. Using a simple 
1/V max (Schmidt 1968; Tinney et al. 1993; Mendez 2000) analysis, the white dwarf number 
density in the direction of the HUDF can be computed. The detection of 4 white dwarfs in 
the HUDF implies a local density of (3.5 ± 1.5) x 10~ 5 ^ stars/pc 3 ^ (1.1 ± 0.3) x 10" 2 . 
The upper and lower limits on the density are computed using old or young white dwarfs 
respectively, while the errors reflect the large uncertainties in the intrinsic luminosities of 
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young and old white dwarfs (~ 1.0 mag). Note that this estimate is actually not affected 
by whether or not object 9230 was included as a HUDF white dwarf. The contribution 
of this object to the 1/V max is negligible. For comparison, Liebert, Dahn & Monet (1988) 
previously determined the local white dwarf density to be 3.0 x 10~ 3 stars/pc 3 , while Reid, 
Sahu & Hawley (2001) found a value of (3.26 ± 1.23) x 10~ 3 stars/pc 3 , and Mendez (2000) 
derived a density of white dwarfs in the thick disk of (2.610 ± 0.59) x 10~ 2 stars/pc 3 . The 
white dwarf density derived here is consistent with most of the previous density estimates. 
Properly determining the disk and halo membership of the 4 white dwarfs identified in the 
HUDF will be needed before a stronger conclusion can be made. Re-observing the HUDF 
field in the i 775 band in 1.5 years would allow (assuming that a 0.18 HUDF pixel shift corre- 
sponds to a 3cr level detection) one to unambiguously detect the proper motion of halo white 
dwarfs with Vr ^ 100 km/s at distances up to 1000 pc. 

The above density estimate assumed that all four HUDF white dwarfs were Galactic disk 
objects. Could some of these objects be in the Galactic halo? At least two of these objects 
(4839, 7768) have Vr ma x values that are lower than 60 km/s. The probability that this is 
caused by a projection effect is under 10% (assuming a relatively low velocity of ~ 100 km/s 
for halo objects). The last two white dwarfs listed in Table 3 (4322, 9020) are at slightly 
larger distances which average to 583 pc, and more importantly have larger values of Vrmax- 
The average distance of these objects still seems a bit low for objects which would be part 
of a Galactic halo whose density increases all the way up to 7 kpc (Reid 1993), while white 
dwarfs with absolute magnitudes of My 606 = 17 mag (Harris et al. 2003) should be detected 
all the way out to 1000 pc. Unless these objects are intrinsically much dimmer, one would 
expect the average distance to these two objects to be larger than 583 pc and much closer to 
our limit of lOOOpc. The Vrmax values inferred for objects 4322 and 9020, while higher than 
those of objects 4839 and 7768, fall substantially short of expected typical halo velocities. 
The probability of observing each of the objects with the V T max values listed in Table 3 is 
between 14-20%. As a group, the probability that objects 4322 and 9020 are halo objects is 
under 3%. Similarly, the hypothesis that all 4 detected white dwarfs in the HUDF are halo 
rather than disk objects is excluded at the 99.9986 % level. Overall, the distances to the four 
white dwarf candidates identified here is more consistent with them being part of a Galactic 
thick disk. Majewski & Seigel (2002) derived a disk scale height of 400-600 pc using a study 
of low velocity white dwarfs. This value was somewhat larger than previous estimates of 
250-350 pc but it is interesting to note that the M-dwarf disk scale height derived in Section 
5.1, as well as the distance estimates of the white dwarfs listed in Table 3 all appear to agree 
with this value. 

It would be interesting to attempt to set an upper limit on the halo white dwarf density based 
on the finding of up to two white dwarfs in the HUDF. Under this assumption, the upper 
limit to the Galactic halo white dwarf density is computed to be (6.53 ±3) x 10~ 3 stars/pc 3 . 
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This result is consistent with the previous work of Mendez & Minniti (2000) using the com- 
bined observations of the Hubble Deep Fields North and South (7.73 x 10~ 3 stars/pc 3 ). 
A population II halo white dwarf population that is older than 12-13 Gyr would have re- 
mained undetected here (An intrinsically fainter than My 606 = 17.0 mag white dwarf pop- 
ulation would also help explaining the low average distances of the white dwarfs listed in 
Table 3). As discussed above, the white dwarfs identified in this work are more likely to be 
within the Galactic thick disk. It is to be noted that Kilic et al. (2004) recently observed 
a disproportionately large fraction of disk to halo white dwarfs in the Hubble Deep Field 
North (HDF-N). Kilic et al. (2004) found no evidence of any white dwarfs with tangential 
velocities larger than 30 km/s down to i 775 = 27.5 mag and the authors concluded that the 
blue objects they saw were all part of the Galactic disk. Similarly, the HUDF observations 
presented here successfully identified 4 white dwarf candidates with i 775 < 27.0 mag in the 
HUDF, none very likely to be in the Galactic halo while 2-3 detection would have been 
expected to be consistent with previous studies. 



5.2.3. Dark halo white dwarfs? 

While dynamical studies of the Galaxy predict the dark matter white dwarf density as 
high as 1.26 x 10~ 2 stars/pc 3 (assuming that 100% of the dark matter halo mass is in the 
form of 0.6M© white dwarfs), the failure to detect a significant population of high velocity 
white dwarfs in the HUDF points to a dark matter halo devoid of a significant white dwarf 
population. Following the methodology of Kawaler (1996), one can estimate the amount 
of dark matter halo mass that is probed by the HUDF images, and directly compute the 
fraction of the dark matter halo mass that could be contained in a white dwarf population. 
Assuming a white dwarf absolute M i775 magnitude of 17.0 mag(Harris et al. 2003) and a 
limiting magnitude of i 77 5 = 27.0 mag, the maximum probed distance is 1000 pc. In the 
corresponding volume, and in the direction of the HUDF, the included dark matter halo 
mass is 3 M Q . A population of dark matter halo white dwarfs could possibly have remained 
undetected in this study for several reasons: (1) The method used in this paper may be 
inefficient at properly identifying halo white dwarfs in Table 3. In this case, if one assumes 
that the HUDF contains 4 halo white dwarfs (i 775 ^ 27.0), these would account for 2M Q , 
a significant fraction of the expected dark matter halo mass. This would however increase 
the ratio of halo to disk white dwarfs to unrealistically high levels (i.e. 3> 2%), and would 
assume that all the white dwarfs seen in the HUDF are from the dark matter halo and 
not from the pop II halo, which is highly unrealistic. The lack of proper motion detections 
also makes this scenario unlikely; (2) White dwarfs do not contribute significantly to the 



- 15 - 



Galactic dark matter halo mass. From microlensing experiements towards the Magellanic 
clouds, the contribution of MACHOs to dark matter has been estimated to be anywhere 
between 20% to less than 2% (Alcock et al. 2000; Afonso et al. 2003; Sahu 2003; Evans 
& Belokurov 2004, and references therein). Furthermore, Majewski & Seigel (2002) showed 
how an increased thick disk scale height of 400-600 pc makes it unlikely that white dwarfs 
could be the MACHO objects; (3) White dwarfs could contribute significantly to the mass 
of the Galactic dark matter halo but are too faint to be detected in this study. Richer et al. 
(2000) showed that an increase in the age of the white dwarf halo population from 10 Gyr 
to 16 Gyr is expected to reduce the number of white dwarf detection (down to V 60 6 = 28.0) 
by a factor of 7. The lack of high velocity white dwarf detections in the HUDF puts upper 
limits on the contribution of an hypothetical white dwarf population to the Galactic dark 
matter halo. Based on the brightest objects in the HUDF for which we have spectra (i.e. 
i 775 ^ 27.0), if the dark matter halo is as old as about 12 Gyr (Hansen 2002; Charboyer et 
al. 1998) then dark matter halo white dwarfs with My = 17.5 (M; 775 = 17.0) contribute less 
than 20% to the dark matter. If the age of the dark matter halo is 10 Gyr, white dwarfs 
(M i775 = 15 . 9) contribute less than 2% to the dark matter. 

The possibility that there is a large population of faint white dwarfs that has remained 
unidentified in this study and which would account for a significant fraction of the dark 
matter halo can be investigated a little further. Recall that in Section 4.3, 14 sources with 
(27.0 ^ i 775 ^ 29.5) and blue ((B 435 — V 60 6) ^ 1-8, (V 60 6 —1775) ^ 0.8) colors were identified. 
Not all of these sources were expected to be bona-fide stars and some are likely to be extra- 
galactic sources, but as discussed in Section 3.1, it is estimated that 95% of the stars in 
the HUDF were properly identified, so that this number should be considered to be an 
upper limit on the number of faint starts in the HUDF. Following the same methodology 
used above, but now using a limiting magnitude of i 77 5 ^ 29.5, the HUDF images probe 
through 79 M Q of the dark matter halo. Even under the extreme assumption that all 14 
unresolved candidates in the UDF are high velocities white dwarfs, this implies that, when 
reaching down to i 775 =29.5, faint white dwarfs in the Galactic dark matter halo contributing 
less than 10% to the total dark matter halo mass for a dark matter halo age of 12 Gyr. 
Future observations of the HUDF would allow more sensitive measurements of the proper 
motion of these faint objects (Section 5.2.2). Setting tighter constraints on the maximum 
tangential velocities of these faint HUDF unresolved objects would allow to exclude some of 
these objects from high velocity white dwarfs and would allow to further lower the maximum 
fraction of the dark matter halo mass that could be accounted for by white dwarfs. 
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6. Conclusions 

A systematic search for unresolved objects in the HUDF identified 46 objects with 
ij75 ^ 29.5 mag. Using the GRAPES spectra of these objects, 28 objects with 1775 $C 27.0 mag 
were spectroscopically identified, including 18 M and later dwarfs (including 2 unconfirmed L 
dwarf candidates), and 2 QSOs. The M dwarf luminosity function of M4 and later-type stars 
was computed and compared to predictions based on Monte-Carlo simulations. Assuming a 
simple analytical model of the Galactic disk, the number of detected M dwarfs was shown 
to be large enough to constraint the scale height of the M and L dwarf Galactic disk to be 
ho = 400 ± lOOpc. Out of 8 stars that were found to not be M or later-type stars, four 
stars were determined to be old white dwarf candidates. Not a single object was found 
to have a proper motion that is larger than 0.027" /year making them likely Galactic thick 
disk objects. Further imaging of the HUDF, using a time span larger than 1.5 year would 
provide more sensitive proper motion measurements. More sensitive proper measurements 
are required to positively place unresolved objects in the Galactic disk or halo. It would also 
allow to search for high tangential velocity objects amongst the 14 fainter (27.0 < i 775 ^ 29.5) 
unresolved sources identified in the HUDF. Excluding halo membership of objects down to 
i775 = 29.5 based on proper motion measurements would further constraint the maximum 
white dwarf contribution to the dark matter halo mass. The currently available observations 
of the HUDF, spanning 73 days, show the absence of a significant population of high velocity 
white dwarfs down to i 775 = 27.0, and a relatively small number of unresolved, faint blue 
sources in the field down to 1775 = 29.5. This is interpreted as a consequence of white dwarfs 
accounting for less than 10% of the dark matter halo, assuming that the dark matter halo 
is 12 Gyr old. 

We would like to thank I. N. Reid for helpful discussions. This work was supported 
by grant GO -09793. 01-A from the Space Telescope Science Institute, which is operated by 
AURA under NASA contract NAS5-26555. 
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Table 1. Unresolved objects in the HUDF, as selected using the criterion <S zg50 ^0.05 

(F850LP ACS band). 
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Table 1 — Continued 
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B435 


B435-V6O6 


V6O6-1775 


1775-Z850 


^1775 


^850 


9006 


53.1853179 


-27.7799976 


27.46 


-0.05 


-0.07 


0.24 


0.039 


0.024 


9020 


53.1685421 


-27.7805214 


27.00 


0.28 


-0.01 


-0.16 


0.025 


0.034 


9212 


53.1485421 


-27.7701387 


25.50 


1.82 


0.96 


0.29 


0.011 


0.003 


9230 


53.1580164 


-27.7691869 


21.18 


0.76 


0.38 


0.10 


0.037 


0.007 


9331 


53.1638484 


-27.767128 


27.95 


2.20 


1.84 


0.67 


0.025 


0.005 


9351 


53.1781506 


-27.7691241 


26.21 


2.03 


1.32 


0.42 


0.015 


0.004 


9397 


53.162852 


-27.7671662 


21.44 


0.26 


0.11 


0.02 


0.043 


0.033 


9959 


53.1611698 


-27.7555109 


29.12 


2.42 


0.91 


0.29 


0.022 


0.021 
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Table 2. Spectral types of the 46 unresolved objects in the HUDF (28 brighter than 
i 775 = 27.0 mag) and their distance estimates assuming that they are main-sequence objects 
or white dwarfs with My = 14 or My = 17. Two additional, fainter extra-galactic sources 
are included (UID 4120, 8157) The spectral types were determined by fitting the GRAPES 
(Pirzkal et al. 2004) spectra of these sources to a set of templates (See Section 4.1). 



UID 


e 

^typc 


i 775 > 27.0 


Mv 606 


MS 


Log(D) 
WD M v = 14 


WD My = 17 


19 


M1-M2 


no 


9.30 ±0.25 


4.072 ±0.050 


3.129 


2.529 


366 


L0-L1 


no 


18.9 ±0.50 


2.780 ±0.100 


3.750 


3.150 


443 


L0-M6 


no 


15.6 ±2.9 


4.261 ±0.570 


4.274 


3.674 


834 


M4-M5 


no 


11.4 ±0.39 


3.777 ±0.077 


3.258 


2.658 


911 


M4-M5 


no 


11.4 ±0.39 


3.338 ±0.077 


2.818 


2.218 


1147 


F6-F8 


no 


3.82 ±0.17 


4.180 ±0.033 


2.142 


1.542 


1343 


A0-K7 


yes 


4.32 ±3.58 


6.080 ±0.716 


3.713 


3.113 


2150 


M3-M4 


no 


11.3 ±0.23 


2.862 ± 0.045 


2.315 


1.715 


2368 


09-M6 


yes 


4.19 ±8.55 


7.145 ±1.710 


3.774 


3.174 


2457 


M0-M1 


no 


8.77 ±0.28 


4.258 ±0.056 


3.209 


2.609 


2977 


A0-M4 


yes 


11.28 ±0.23 


4.365 ±0.045 


3.818 


3.218 


3166 


K4-K5 


no 


6.96 ±0.14 


3.573 ±0.027 


2.164 


1.564 


3561 


M0-M1 


no 


8.77 ±0.28 


4.697 ±0.056 


3.647 


3.047 


3794 


M0-M1 


no 


8.77 ±0.28 


4.594 ±0.056 


3.544 


2.944 


3940 


05-M3 


yes 


2.37 ±7.94 


7.315 ±1.588 


3.702 


3.102 


4120 


GAL 


yes 










4322 


A7-K5 


no 


4.70 ±2.4 


5.632 ± 0.479 


3.548 


2.948 


4643 


A0-M6 


yes 


6.75 ±5.96 


6.578 ± 1.191 


4.237 


3.637 


4839 


G5-K2 


no 


5.65 ±0.61 


5.087 ±0.122 


3.400 


2.800 


4945 


M4-M5 


no 


11.4 ±0.38 


3.315 ±0.077 


2.796 


2.196 


5317 


A0-M6 


yes 


6.71 ±6.0 


6.267 ± 1.199 


3.910 


3.310 


5441 


M2-M3 


no 


9.93 ±0.38 


4.228 ±0.076 


3.407 


2.807 


5921 


K4-K5 


no 


6.96 ±0.14 


3.565 ±0.027 


2.155 


1.555 


5992 


K7-M0 


no 


8.07 ±0.43 


4.568 ±0.085 


3.373 


2.773 


6334 


F8-L0 


yes 


8.34 ±4.4 


5.505 ±0.875 


3.791 


3.191 


6442 


A0-G5 


yes 


2.89 ±2.2 


6.271 ±0.429 


3.865 


3.265 


6461 


M0-M1 


no 


8.77 ±0.28 


4.251 ±0.056 


3.202 


2.602 


6620 


B8-M4 


yes 


5.45 ±5.6 


6.345 ±1.121 


3.812 


3.212 


6732 


GAL 


no 










7113 


G8-K0 


yes 


5.60 ±0.19 


5.543 ±0.038 


3.861 


3.261 


7194 


05-K5 


yes 


0.76 ±6.3 


7.215 ±1.266 


3.601 


3.001 


7357 


09-M5 


yes 


3.73 ±8.1 


7.131 ± 1.618 


3.760 


3.160 


7525 


M1-M2 


no 


9.33 ±0.23 


4.429 ± 0.045 


3.491 


2.892 


7768 


K5-K7 


no 


7.37 ±0.27 


4.616 ±0.055 


3.286 


2.686 


7894 


B8-K7 


yes 


3.73 ±3.9 


6.223 ±0.780 


3.680 


3.080 


8081 


05- A3 


yes 


-1.98 ±3.6 


7.471 ±0.716 


3.843 


3.243 


8157 


GAL 


yes 
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Table 2 — Continued 



UID 


Stype 


i 7 75 > 27.0 


M v 606 




Log(D) 










MS 


WD My = 14 WD My = 17 



8186 


B0-K7 


yes 


1.89 ±5.8 


7.038 ± 1.151 


3.764 


3.164 


9006 


B9-A3 


yes 


0.94 ±0.65 


6.328 ±0.131 


3.697 


3.097 


9020 


F0-F2 


no 


3.16 ±0.44 


5.716 ±0.089 


3.540 


2.940 


9212 


M0-M1 


no 


8.77 ±0.28 


3.979 ±0.056 


2.929 


2.329 


9230 


K2-K3 


no 


6.40 ±0.14 


3.800 ±0.028 


2.279 


1.679 


9331 


M5-M6 


no 


12.3 ±0.45 


3.701 ± 0.089 


3.345 


2.745 


9351 


M1-M2 


no 


9.30 ±0.25 


3.974 ±0.050 


3.031 


2.431 


9397 


GAL 


no 










9959 


M0-M1 


no 


8.77 ±0.28 


4.585 ±0.056 


3.535 


2.935 



Table 3. White dwarf candidates in the HUDF (i 775 ^27.0 mag). We list the derived 
distance to each object under the assumptions that it is a main sequence star or white 

dwarfs (Richer et al. 2000) with the appropriate V — % colors (See Table 1). Vrmax is the 
upper limit of the tangential velocity (Vr) of each object. The seventh column lists the 
acceptable nature of the object while the last column gives the l/V max density for that 

particular object /scenario. Section 5.2 explains how these quantities were derived in more 

detail. 



UID 


Type 


v 606 


D 


Age 


Vr 


Possible 


/ * 1 1 L LI Jl, 










(pc) 


(Gyr) 


(km/s) 


Type 


(star /pc~ 


3 ) 

/ 


1147 


MS F6-F8 


3.7-4.0 


1.5 ± 1.1 x 10 4 






F6-F8 






1147 


Young WD 


14.0-15.0 


110 ± 24 


1.4-3.9 


14.1 








1147 


Old WD 


17.4-18.3 


24 ± 5 


10.3-12.1 


3.0 








3166 


MS K4-K5 


6.8-7.1 


3.7 ±0.23 x 10 3 






K4-K5 






3166 


Young WD 


14.8-15.8 


80± 18 


2.5-5.7 


10.2 








3166 


Old WD 


17.2-18.2 


27 ± 6 


9.9-11.5 


3.5 








4322 


MS A7-K5 


2.3-7.1 


4.29 ±3.4 x 10 5 












4322 


Young WD 


12.1-13.1 


6.9 ± 1.6 x 10 3 


0.4-1.4 


1545.6 








4322 


Old WD 


17.5-18.4 


590 ± 126 


10.7-12.9 


75.5 


Disk or Halo 


3.37 x 10 


-3 


4839 


MS G5-K2 


5.0-6.3 


1.2 ±0.3 x 10 5 












4839 


Young WD 


12.5-13.5 


2.4 ±0.5 x 10 3 


1.1-3.3 


885. 


Disk 


9.45 x 10 


-6 


4839 


Old WD 


17.4-18.3 


417 ± 88 


10.4-12.3 


53.4 


Disk or Halo 


1.80 x 10 


-2 


5921 


MS K4-K5 


6.8-7.1 


3.7 ±0.2 x 10 3 






K4-K5 






5921 


Young WD 


14.5-15.5 


93 ±20 


1.9-4.7 


11.9 








5921 


Old WD 


17.3-18.2 


26 ±5 


10.1-11.7 


3.3 








7768 


MS K5-K7 


7.1-7.6 


4.1 ±0.5 x 10 4 












7768 


Young WD 


14.5-15.5 


1.3 ±0.3 x 10 3 


1.9-4.7 


160.2 


Disk 


2.1 x 10- 


-5 


7768 


Old WD 


17.3-18.2 


345 ± 70 


10.1-11.7 


44.1 


Disk or Halo 


9.9 x lO - 


-4 


9020 


MS F0-F2 


2.7-3.6 


5.2 ± 1.0 x 10 5 












9020 


Young WD 


13.3-13.3 


6.2 ± 1.4 x 10 3 


0.5-1.6 


793.9 








9020 


Old WD 


17.5-18.4 


576 ± 123 


10.7-12.9 


73.7 


Disk or Halo 


3.16 x 10 


-3 


9230 


MS K2-K3 


6.3-6.7 


6.3 ±0.4 x 10 3 






K2-K3 






9230 


Young WD 


13.6-14.6 


187 ± 41 


1.1-3.2 


23.9 








9230 


Old WD 


17.5-18.4 


32 ± 7 


10.4-12.3 


4.0 
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Magnitude 



Fig. 1. — The distribution of <5>i 775 (green) and <S Z850 (red) values for the 10181 HUDF objects. 
The stellarity index S was defined (see Section 3) so that higher values correspond to objects 
which are increasingly resolved in the images. We selected unresolved objects fainter than 
i 775 ^ 20.0 mag to be unresolved if 5 zg50 ^0.05, and <S Z850 ^0.15 for brighter objects which 
are saturated. 
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Fig. 2. — The distribution of <S Z850 values for the brightest 369 objects in the field. Each 
object appears 10 times in this figure, with increasingly lower signal to noise (increasingly 
darker dots). The observed increase in S values as a function of magnitude at i 775 ^ 29.0 mag 
is similar to the one observed in Figure 1 where the real distribution of the S values from the 
10181 HUDF objects is shown. This demonstrates the robustness of S to distinguish between 
point sources and extended objects down to faint magnitudes. At fainter magnitudes, images 
of stars become increasingly dominated by noise which tends to increases the measured values 
of S . The same test was performed in the i 7 7 5 band and yielded an identical behavior. 
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Fig. 3. — GRAPES spectra of unresolved objects in the HUDF and the two best matching 
templates (solid lines) listed in Table 2 
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Fig. 4. — GRAPES spectra of unresolved objects in the HUDF and the two best matching 
templates (solid lines) listed in Table 2 



-29- 




6000 7000 8000 9000 6000 7000 8000 9000 

Wavelength A 



Fig. 5. — GRAPES spectra of unresolved objects in the HUDF and the two best matching 
templates (solid lines) listed in Table 2. The bottom 4 objects are the extra-galactic sources 
4120, 6732, 8157, 9397 estimated to be at redshifts of z = 2.1, z = 3.2, z ^ 3.0, z = 3.0, 
respectively (Xu 2004). 
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Fig. 6. — The location of the 28 spectroscopically identified sources in a magnitude color plot. 
The small, dark points are the entire set of objects in the HUDF catalog with i 775 29.0 
mag and z 850 ^ 28.2 mag. Solid red squares are M dwarfs. Empty red squares are the two L 
dwarf candidates. Solid green triangles represent the non M or L dwarfs stars which could 
be main-sequence stars. The empty blue circles are white dwarf candidates (Table 3). The 
solid yellow circles represent the QSOs. 
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7. — The location of the 28 spectroscopically identified sources in a magnitude color 
The symbols are the same as in Figure 6. 
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Fig. 8. — Left panel: observed shifts (scaled by a factor of 500) of the 360 most compact 
sources in the HUDF between Epoch 1 to Epoch 2, i 775 band, images; Right panel: cor- 
rected shifts, obtained after applying a correction based on a fit of the raw measurements 
(right panel, see Section 3.2). North is up on these panels. The axes units are UDF pixels 
(0.030" /pixel). 
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Fig. 9. — (V6O6-1775) vs - (B435-V6O6) color-color plot of 1775 ^ 27.0, unresolved, spectroscop- 
ically identified objects in the HUDF. M and L dwarfs are shown using red squares with 
photometric error bars. Main-sequence stars are shown in green thick circles. Extra-galactic 
objects, including those dimmer than 1775 = 27.0 mag, are shown using large blue stars. The 
5 white dwarf candidates are shown using large blue crosses. Pickles main-sequence objects 
are also shown using small black circles. The (B 435 -V606) locations of the O, A , F, G, K, and 
M stellar type are shown for reference. Hot white dwarfs are shown in large blue squares. 
Cool white dwarfs ((3000°K, 4000°K, 5000°K), Harris et al. 2001) are shown using labeled, 
empty circles. Finally, the solid black line shows the QSO track, from redshift of z=4.0 (top 
right) to z=0.1 (bottom left) 
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Fig. 10. — (1775-Z850) vs - (^606-1775) color-color plot of i 775 ^ 27.0 mag, unresolved, spectro- 
scopically identified objects in the HUDF. The labels and symbols used are the same as in 
Figure 9. 
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Fig. 11. — Predicted and observed cumulative luminosity functions of stars with spectral 
types M4 and later. This plot shows the expected distribution from the Galactic disk and 
halo. Four different disk scale heights are plotted with h Q values of 200, 300, 400 and 
500 pc (bottom curves to top curves respectively). For each value of ho, Mass Functions 
corresponding to a = 0.5 and 1.5 are shown (solid and dashed lines respectively). The 
observed HUDF distribution is shown using a solid black line which is best fit by a disk with 
a scale height of ho = 400pc while reasonably excluding values of h = 200,300, and 500. 



